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I. INTRODUCTION 


Applications of perturbative QCD to deeply virtual Compton scattering (|l|-Q and hard exclusive electropro¬ 
duction processes ili require a generalization of usual parton distributions for the case when long-distance 
information is accumulated in nonforward matrix elements {p — r\O{0, z)\p)\z 2 ^Q of quark and gluon light-cone 
operators. In refs. it was shown that such matrix elements can be parametrized by two basic types of non- 

perturbative functions. The double distribution F{x,y;t) specifies the light-cone “plus” fractions xp'^ and ?/r+ of 
the initial hadron momentum p and the momentum transfer r carried by the initial parton. Since r'^ is proportional 
to p'^\ r+ = , it is possible to introduce the nonforward parton distribution Tq {X ; t) with X = x + yC, being 

the total fraction of the initial hadron momentum taken by the initial parton^. For processes mentioned above, 
the parameter C = 1 — {p'z)l{pz) characterizing the longitudinal momentum asymmetry (“skewedness”) of the 
nonforward matrix element takes the values 0 < ^ < 1- 

At leading twist, there are two light-ray quark operators '(/i(0)7^if(0, z; and '0(0)7^75if(0, z; 

where E{0,z;A) is the standard path-ordered exponential which makes the operators gauge-invariant. In the 
forward case, the first operator is related to the spin-averaged distribution functions fa{x) while the second one 
corresponds to the spin-dependent distribution functions Afa{x). The nonforward parton distributions related 
to the ip{0)j^E{0, z; A)'tp{z) operators were studied in refs. In this paper, we will discuss flavor-singlet 

parton helicity-dependent nonforward distributions corresponding to quark operators ' 0 ( 0 ) 7 ^ 75111 ( 0 , z; A)'0(z) and 
the gluonic operator G)j,a{0)E{0, z', A)Gaiy{z) mixing with each other under evolution. 


II. NONFORWARD DISTRIBUTIONS 

We define the nonforward quark distributions by writing the relevant matrix element as (cf. [|^,Q) 

(p',s'|0a(O)z75£l(O,z;A)0a(2) b,s)U2^o (1) 

= u{p\s')zj5u{p,s) + dX 

+ M u{p', s'h,u{P, s) I" t) + 1)) dX, 

where t = {p' — p)^, a denotes the quark flavor (here we consider only the flavor-diagonal distributions), M is 
the nucleon mass and s, s' specify the nucleon polarization. Throughout the paper, we use the “hat” convention 
z = z^ 7 ^. InEq. i) , we explicitly separated quark and antiquark contributions (cf. [Q). This definition corresponds 
to the parton picture in which the initial quark (or antiquark) takes the momentum Xp from the hadronic matrix 
element and “returns” into it the momentum {X — ()p. Since the fraction X — C is positive for X > C and negative 
when X < 0, the nonforward distributions can be divided into two components. In the region X > 0, one can treat 
fy“(X, t) as a generalization of the usual distribution function Afa{x). In particular, in the limit t —> 0, 0 —> 0, the 
limiting curves for Qq{XA) reproduce A/a(X): 

g^^oiX,t = 0) = Afa{X) ; ^“^o(X,t = 0)=A/a(X). (2) 

On the other hand, in the region X < 0, both quarks should be treated as going out of the nucleon matrix element, 
with momenta Xp and (0 — X)p, respectively. Now, one can define Y = X/0 and treat the function t/“(X) as the 
distribution amplitude 'I'“(X). In particular, the ^-part in this region can be written as 


*The off-forward parton distributions introduced by X. Ji |i|Q (see also ^) and non-diagonal distributions of Collins, 
Frankfurt and Strikman [Q can be related to nonforward distributions (see ^ ) but do not coincide with them. 
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Cu{p')zu{p) [ + dr = Cu(y)^M(p) / (3) 

Jo ^ Jo 

where the distribution amplitude d>^(F) is defined by = Q^{Y() + Q^{Y(). The function gives the 

probability amplitude that the initial nucleon with momentum p is composed of the final nucleon with momentum 
p' = p — r and a qq pair in which the pair momentum r is shared in fractions Y and 1 — Y = Y. 

For gluons, the nonforward distribution is defined through the matrix element 


{p' I z^z.G;^i0)E^\0, z; A)GUz) I p)U-o 
/■i 


( 4 ) 


= u{p')zj5u{p) (z-p) J 2 L 

+ U{p')^j5uip){z -p) ^ 


g-iX(pz) _ ^i{X-Q(pz) 


2 L 


^-iX(pz) _ gi(X-C)(pz) 


gl{x-t)dx 

Vl{X-t)dX. 


, ^av — 2 ^avi 3 pG^^. Since there are no “antigluons”, the exponentials e ^x{pz) gz(x C.){pA gj.g 


As usual, Gai, = 

accompanied here by the same function {X\t). Again, the contribution from the region 0 < AT < C can be 


written as 


iu{p')z-f^u{p) (z-r) [ t) dY + “V” term, 

Jo 

with the Y <-^Y antisymmetric generalized distribution amplitude given by 


( 5 ) 


( 6 ) 


In the formal t = 0 limit, the nonforward distributions g^(X;t), 'P^{X]t) convert into the asymmetric distribution 
functions g^{X), V^{X). Finally, in the C = 0 limit, G^iX) reduces to the usual polarized gluon density 


GU^X) = XAg{X). 


Under pQCD evolution, the gluonic operator 


Og{uz,vz) = Zg,z^G'ta{uz)E°‘^(uz,vz] A)&^^{vz) 


mixes with the flavor-singlet quark operator 


N; 


Oq{uz, vz) = O^'^^uz, vz) 


a=l 


where 


0^+\uz,vz) = i 


lj)a{uz)z-ioE{uZ, VZ] A)'lj)a{vz) + tpaivz)z"f5E{vZ, UZ] A)tjja{uz) 


( 7 ) 


( 8 ) 


(9) 


( 10 ) 


The nonforward distribution function g9 {X ; t) for the flavor-singlet quark combination 


(p',s'| Oq{uz,vz) b,s)U2^o = u{p’, s')zj 5 u{p, s) 


2 L 


^—ivX{pz)-\-iuX\pz) _j_ ^ivX' {pz)—iuX{pz) 


(where X' = X — Q can be expressed as the sum of “a H- a” distributions: 


Nf 




G^ {X]t) dX + “P” term, 

( 11 ) 


( 12 ) 
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Writing the contribution from the 0 < X < C region as 


Cu{p')z-r 5 u{p) (z-r) f t) dY + “F’ term, (13) 

Jo 

we introduce the flavor-singlet quark distribution amplitude (Y ; t) which has the symmetry property = 

(Y ; t) with respect to the Y ^ Y transformation. 


III. EVOLUTION EQUATIONS FOR LIGHT-RAY OPERATORS 


Near the light cone ~ 0, the bilocal operators 0{uz,vz) develop logarithmic singularities Inz^. Calculation- 
ally, these singularities manifest themselves as ultraviolet divergences for operators taken on the light cone. The 
divergences are removed by a subtraction prescription characterized by some scale fi: Q^[X]t) p). At 

one loop, the set of evolution equations for the flavor-singlet light-ray operators has the following form (cf. [p|,p^): 

d 

p — Oa{0,z)= Aab{u,v)Ob{uz,vz)6{u +V < 1) dudv , 

Jo Jo V 


(14) 


where a,b = g,Q and v = l— v, u = l — u. For flavor-nonsinglet distributions, there is no mixing, and their 
evolution is generated by the QQ-kernel alone. To calculate the kernels, we incorporated the approach JK] based 
on the background-field method. Below we present our results in the form similar to that used in refs. |6 


^QQiux) = — Cf ( 1 + -S{u)S{v) -f S{u) 




-I- {u •*-> u} 


AgQiu, ^ Cf ^(5(m)(5(z;) - 2^ , 
v) = —Nf{l-u-v) , 

TT 


Agy(u,v} = '^N^[ 4{1 - u - v) + X S{u}S{v) + -J 6{u) 


^ - S(.) t A 


V 


'0 


V 


+ {u tj} 


(15) 

(16) 

(17) 

(18) 


Independently, these kernels were calculated by Blumlein, Geyer and Robaschik ||^,|^. Their results agree with 
ours. 


IV. EVOLUTION EQUATIONS FOR NONFORWARD DISTRIBUTIONS 

Inserting the light-ray evolution equations ( p^ between chosen hadronic states and parametrizing matrix ele¬ 
ments by appropriate distributions, one can get the “old” DGLAP nn and BL-type @0 evolution kernels 
as well as calculate the new kernels AIV“^(A, Z) governing the evolution of nonforward parton distributions: 

p±g-{X;t-,p)= [ y^AW^\X,Z)gyZ;t-,p)dZ. (19) 

Jo ^ 

Extracting AW^^{X, Z) from the light-ray kernels Aab{u,v), one should take into account the extra (pz) factor 
in the rhs of the gluon distribution definition, which under the Fourier transformation with respect to (pz) results 
in the differentiation djdX. Thus, it is convenient to introduce first the auxiliary kernels AM“*’(A, Z) directly 
related to the light-ray kernels A{u,v) by 

AMq^{X,Z)= f [ Aab{u,v) S{X — uZ + v{Z — ()) 6{u + V < 1) dudv . (20) 

Jo Jo 
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The Ahh-kernels are obtained from the AM-kernels using 

AtT®®(A, Z) = AM^^X, Z) , AW^^{X, Z) = AM^^{X, Z), (21) 

Z) = -AMf^{X, Z) dX , AW^^X, Z) = AM^®(A, Z). (22) 

Hence, to get AH(?^(A, Z) we should integrate AM®‘^(A, Z) with respect to X. We fix the integration constant 
by the requirement that AW^*^(A, Z) vanishes for A > 1. Then 


AW^^{X,Z) = J AMf^{X,Z)dX. 


(23) 


Integrating the delta-function in eq.(|20|) produces four different types of the 0-functions, each of which corresponds 
to a specific component of the kernel governing the evolution of the nonforward distributions. 


V. BL-TYPE EVOLUTION KERNELS 

When C, = 1, Gc{X) reduces to a distribution amplitude whose evolution is governed by the BL-type kernels: 

A (A, Z) = 1/“'' (A, Z) . (24) 

Taking ^ = 1 in Eq.(^) we obtain 

AM“^i(A,Z) = t/“^(A,Z) = / f Aab{u,v)S{X -uZ-v{l-Z))9{u + v <l)dudv. (25) 

^0 Jo 

In fact, the BL-type kernels appear as a part of the nonforward kernel IU“^(A, Z) even in the general ( ^ 1,0 
case. As explained earlier, if A is in the region A < C, then the function Gc{X) can be treated as a distribution 
amplitude 'I',^(Y) with Y = A/^. For this reason, when both A and Z are smaller than C, the kernels H)?^(A, Z) 
simply reduce to the BL-type evolution kernels V°'^{X/C,,Z/Q. Indeed, the relation ( ^ ) can be written as 

AMf{X,Z) = j [ [ Aab{u,v)S{X/(:-uZ/(;-v{l-Z/(:)) e{u + v<l)dudv. (26) 

C Jo Jo 

Comparing this expression with the representation for the {7“^(A, Z) kernels, we conclude that in the region where 
A/C < I and Z/( < 1, the kernels AM“*’(A, Z) are given by 

AMf (A, Z)\o<[x,z}<c = ^ {X/C, Z/C) ■ (27) 

Now, using the expressions connecting the AW- and AM-kernels, we obtain the following relations between the 
nonforward evolution kernels AM//^(A, Z) in the region 0 < {A, Z} < C, and the BL-type kernels U“*'(A, Z)\ 

AW^^{X, ^) = 7 (A/C, Z/Q) ■ AW^^{X, Z) = (A/C, Z/Q ; 

c 

AW^^ (A, z) = 1 yQ® (A/C, Z/0 ; AWf (A, Z) = ^ (A/C, Z/C) . (28) 

The kernels V°'^{X, Z), in their turn, are derived from the auxiliary kernels U°'^{X, Z). Due to the symmetry 
property Aab{u,v) = Aab{v,u) the kernels U°'^{X,Z) satisfy U°‘^{X,Z) = U°‘^{X,Z). Hence, it is sufficient to 
know the [/-kernels in the X < Z region only: 

[/“*'(A, z) = 6»(A < z) u/f'ix, z) + e{z < x)u/;\x, z ), 
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with the basic function Uq^{X, Z) = 9{X < Z) U°^^{X, Z) given by 


U^\X, Z) = ^J^ Aab iv-{X- v)/Z, v) dv . 


Using Eqs.(15)-(|18D, the A ^ Uq conversion formulas 


and Eqs.(p 


6{u) 5{v) ^ S{Z - X) , , 5(m)-^0 , S(u) (-) 

Z V \v/ 


5{v) 


X 


1 


(5(r 


Z ) z-x 

we get the BL-type kernels 




Z Z-X 


u + v 


X 


1 - 


X 


Y^^{x,z) = —Cf 




1 


Z-x 


e{x <z) 


+ {X^X,Z^Z} 


-I + 


yQ9(X,Z) = ^X/ |-J0(X<Z) +£0(X>Z)| 
V^^{X,Z) = ^Cf |^0(X<Z) - ^0(X>Z)|, 


USS(X,Z) = — Xc 

TT 


2X2 _ X - Z 
Z2 


(X < Z) + 


e{x <z) 
z-x 


+ {X ^ X,Z ^ Z} 


-I + 




( 29 ) 


(30) 

(31) 

(32) 

(33) 

(34) 


calculated originally in for flavor-singlet pseudoscalar meson distribution amplitudes. With respect to 

integration over 0 < X <1, the “plus”-prescription for a function U(X, Z) is defined by 


[U(X,Z)]+ = U(X,Z) 


(5 (X - Z) 



u(y, z) dY. 


(35) 


The BL-type kernels also govern the evolution in the region corresponding to transitions from a fraction Z which 
is larger than ^ to a fraction X which is smaller than Y Indeed, using the (5-function to calculate the integral over 
u, we get 


1 ( \ 

AMf(X,Z)|x<c<z = ^y^ [l-X/Z-t-(l-C/Z)],ujdu, (36) 

which has the same analytic form ( |^ ) as the expression for M“*’(X, Z) in the region X < Z < Y For QQ, gg and Qg 
kernels, this automatically means that AIU“^(X, Z)\x<q<z is given by the same analytic expression as A1U“^(X, Z) 
for X < Z < (■. Because of integration in Eq.(p3|), to get AV1(?‘^(X, Z) one should also know AM®*^(X, Z) in 
the region Q < X < Z. However, our explicit calculation confirms that AVE®‘^(X, Z) in the transition region 
X < ((■ < Z is given by the same expression as A1U®'^(X, Z) for X < Z < C,. 

In application to parton distributions related to nonforward matrix elements, X. Ji was the first who calculated 
analogous kernels P'{x,Y) which govern the evolution of his off-forward parton distributions H{x,t;fi) in the 
—^/2 < X < ^/2 region (in our variables this region corresponds to 0 < X < ^). He used a direct momentum- 
representation approach in the light-cone gauge. After proper redefinitions (discussed in [Q), we reproduced his 
expressions for the first three kernels. For the gluon-gluon kernel, our result formally differs from that obtained by 
X. Ji 1^. However, due to the symmetry properties of the gluon distribution in the X. Ji approach, the relevant 
integral vanishes and the difference does not contribute to the evolution. Blumlein et al. derive the “extended” 
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BL-kernels [|| from the light-ray evolution equations. For X ^ Z, we agree with their results except for the 
^Q'kernel and up to obvious misprints in the QQ and g^-kernels |^. 


VI. GENERALIZED DGLAP KERNELS 


When A > ^, we can treat the asymmetric distribution function Q^(X) as a generalization of the usual distri¬ 
bution function A/a (A) for a skewed kinematics. Hence, evolution in the region C<A<l,/<Z<lis close 
to that generated by the DGLAP equation. In particular, it has the basic property that the evolved fraction A 
cannot be larger than the original fraction Z. The relevant kernels are given by 

AMf(A,Z)|c<x<z<i = Aab (w (1 - X/Z), w (1 - X'/Z')) dw, (37) 


where A' = A — / and Z' = Z — C, are the “returning” partners of the original fractions A, Z. Note, that since 
Z — X = Z' —A', the kernels AM^*'(A, Z) are given by functions symmetric with respect to the interchange of A, Z 
with A', Z'. Using the table for transition from the Aab-kernels to the AM“*’-kernels in the region C, < X < Z < \ 


(6{u)-+5{v)- 

\ V u 


S{u) d{v) ^ S{Z - X) ; 1 
1 


Z-A 


A A' 

'z^~x 


Z-X 
ZZ' 


; (u -b u) 


; b{u) -h(5(u) — 

V u 


Z-A 
2ZZ' 

1 


A_^' 
Z Z' 


Z-A 




f- 

\Z' 


and Eqs.(^, (Hj), we obtain the kernels AP/^(A, Z) = AH/‘^(A, Z)|,^<x<z<i: 


ApW(A,Z) = —Cf 

^ TT 


1 


1-b 


AA' 


Z-A 
+ 6{X-Z) 


ZZ' 
3 _ 
2 ” 


6» (A < Z) 


du 

u 


dv 

V 


AP^HX.Z) = + ^AP^,(X/Z ), 


APf{X,Z) = ^CF\^ + ^- 


X X' XX' 


AP®®(A, Z) = 


■Nr_ 


IF 

'X X'' 

\V 

)z (a 


+ Six - Z) 


ZZ' 

Z-X 
ZZ' 

Po 

2Nr. ~ 


|ap,q(A/Z), 


1 


'A 

Z-A I VZ 
du dv 

lo u Jo V 


-APQQiXjZ), 


X' 


0(A < Z) 


^APAXIZ). 


(38) 


(39) 

(40) 

(41) 


(42) 


The formally divergent integrals over u and v provide here the usual “plus”-type regularization of the 1/(Z — A) 
singularities. The prescription following from Eqs.(^),(^^ is that combining the 1/(Z — A) and (5(Z — A) terms 
into [Qq{Z) — t/^(A)]/(Z — A) in the convolution of AP(^{X, Z) with Qq{Z) one should change u —> 1 — A/Z and 
V l-X' jZ'. 

As expected, the AP“*'(A, Z) kernels have a symmetric form. The arrows indicate how the nonforward kernels 
AP“*'(A, Z) are related to the DGLAP kernels in the C = 0 limit when Z = Z' and A = A'. Deriving these 
relations, one should take into account that the gluonic asymmetric distribution function t/®(A) reduces in the 
/ ^ 0 limit to XAg{X) rather than to Ag{X). 


^We are grateful to J. Blumlein who informed us that the authors of ref. Illjl agree with our results. 
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After the appropriate redefinitions, we managed to reproduce from our results all four kernels APab{x, (relevant 
to the X > ^j2 region) calculated by X. Ji Q. 

Note, that in the region Z > ( the evolved fraction X is always smaller than Z. Furthermore, if Z < ^ then 
also X < C, i.e., distributions in the X > C, regions are not affected by the distributions in the X < ^ regions. 
Hence, information about the initial distribution in the Z > C, region is sufficient for calculating its evolution in 
this region. This situation may be contrasted with the evolution of distributions in the Z < C, regions: in that case 
one should know the nonforward parton distributions in the whole domain Q < Z < 1. 


VII. CONCLUSIONS 

In this letter, we discussed the calculation of the evolution kernels AW^{X, Z) for nonforward parton distributions 
Q(^{X,t) sensitive to parton helicities. We presented the evolution kernels for the relevant light-ray operators and 
demonstrated how one can obtain from them the components of the nonforward kernels A.Wc^{X,Z). Our results 
have a transparent relation with DGLAP and BL-type kernels and a compact form convenient for further practical 
applications such as numerical studies of the evolution of nonforward distributions. 
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